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ABSTRACT 

We use the Wide Field Camera 3 onboard the Hubble Space Telescope to obtain deep, 
high-resolution photometry of the young (age ^ 300 Myr) star cluster NGC 1856 
in the Large Magellanic Cloud. We compare the observed colour-magnitude diagram 
(CMD), after having applied a correction for differential reddening, with Monte Carlo 
simulations of simple stellar populations (SSPs) of various ages. We find that the main 
sequence turn-off (MSTO) region is wider than that derived from the simulation of a 
single SSP. Using constraints based on the distribution of stars in the MSTO region 
and the red clump, we find that the CMD is best reproduced using a combination of 
two different SSPs with ages separated by 80 Myr (0.30 and 0.38 Gyr, respectively). 
However, we can not formally exclude that the width of the MSTO could be due to a 
range of stellar rotation velocities if the efficiency of rotational mixing is higher than 
typically assumed. Using a King-model fit to the surface number density profile in 
conjunction with dynamical evolution models, we determine the evolution of cluster 
mass and escape velocity from an age of 10 Myr to the present age, taking into account 
the possible effects of primordial mass segregation. We find that the cluster has an 
escape velocity 14sc — 17 kms“^ at an age of 10 Myr, and it remains high enough dur¬ 
ing a period of ~ 100 Myr to retain material ejected by slow winds of first-generation 
stars. Our results are consiste nt with the presence o f an a ge spread in NGC 1856, in 
contradiction to the results of iBastian fc Silva-Villal (|2013[1 . 
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1 INTRODUCTION 

Over the last ~ half a dozen years, deep colour-magitude 
diagrams (CMDs) from images taken with the Advanced 
Camera for Surveys (ACS) and the Wide Field Camera 3 
(WFC3) aboard the Hubble Space Telescope (HST) revealed 
that several intermediate-age (~ 1 - 2 Gyr old) star clusters 
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A popular interpretation of these extended MSTOs 
(hereafter eMSTOs) is that they 


formed at different times within the parent cluster, with age 
sprea d s of 150-500 Mvr iMlone et al.ll2009l: ICirardi et 5.1 
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2 OI 4 II . Alternative potential causes to explain the eM- 
STO phenomenon presented in the recent literature in- 
dude spreads in rotation velocity amon g turn-off stars 
dBastiaii k, de Mini3l2009l : lYang et al.ll2013l . but see Girardi 
et al. 2011), a pho tometric feature of interact ing binary star s 
(lYang et al.ll20lil') or a combination of both (iLi et al.ll2012l l. 

An important aspect of the nature of the eMSTO phe- 


are due to stars that is not shared by all such clusters (e.g 
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that the key factor in determinining whether or not a cluster 
features an eMSTO is its ability to retain material ejected by 
first-generation stars that feature relatively slow stellar out¬ 
flows (the so-called “polluters”). Following the arguments 
presented in Glia, eMSTOs can only be hosted by clus- 
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ters for which the escape velocity was higher than the wind 
velocity of such polluter stars at the time such stars were 
present in the cluster. 

Currently, the most popular candidates for first- 
generation “polluters” are (i) intermediate-mass AGB 
stars (4 < M/Mp, < 8, hereafter IM-AGB; e.g., 

iD’Antona fc Ventural l2007l . and references therein), (ii) 
rapi dly rotating massive st ars (often referred to as “FRMS”; 
e.g. Decressin_et_^ I [2^ and (Hi) massive binary stars 
dde Mink et al.ll2009l ~). 


For what concerns instead the formation scenario, the 
currently two favored ones predict that the stars have been 
formed from or polluted by gas that is a mixture of pris¬ 
tine material and material shed by such “polluters”. In par- 
ticular, in the “in s i tu star formation” scenario ( see e.g., 
iD’Ercole et al.|[2008l . l201(il : IConrov fc Spergelll201lfl . subse¬ 
quent generations of stars are formed out of gas clouds that 
were polluted by winds of first generation stars to vary¬ 
ing extents, during a period spanning up to a few hun¬ 
dreds of Myr, depending on the nature of the “polluters”. 
Conve rsely, in the altern ative “early-disc accretion” sce¬ 
nario (iBastian et al.l [201,81 1. the chemically enriched mate¬ 
rial, ejected from interacting high-mass binary systems or 
FRMS stars, is accreted onto the circumstellar discs of low- 
mass pre-main sequence (PMS) stars that were formed in 
the same generation during the first ~ 20 Myr after the 
formation of the star cluster. 

Two recent studies p rovided support to the p redictions 
of the “in situ” scenario. iGoudfrooii et all (|2014| . hereafter 
G14) studied in detail HST photometry of a sample of 
18 massive intermediate-age star clusters in the Magellanic 
Glouds that covered a variety of masses, ages, and radii. G14 
found that all the star clusters in their sample host eMSTOs, 
featuring age spreads that correlate with the clusters’ escape 
velocity 14sc, bot h currently and at an age of 10 Myr. 

Furthermore, [^rrenti et al.l (l2014l . hereafter C14) stud¬ 
ied 4 low-mass (« 10"^ M©) intermediate-age star clusters in 
the Large Magellanic Gloud (LMC) and found that the two 
clusters that host eMSTOs have 14sc ^ 15 kms“^ out to 
ages of ~ 100 Myr, whereas 14sc < 12 kms“^ at all ages for 
the two clusters that do not exhibit eMSTOs. These results 
suggest that the “critical” escape velocity for a star cluster 
to able to retain the material ejected by the slow winds of 
the first generation polluters seems to be in the approximate 
range of 12- 15 kms“^. Interestingly, these escape velocities 
are consistent with wind velociti es of IM-AGB stars which 
are in the ra nge 1 2-18 kms“^ (IVassiliadis fc Woodlll993l : 
IZiilstra et ahlllO^ j and with the low end of observed veloc¬ 
ities of ejecta of the massive bi nary star RY Scuti (15-50 
kms~^. ISmith et al.ll200l2007l l. 


However, in the context of the multiple population sce¬ 
nario, one piece of the puzzle that is still missing is the iden¬ 
tification of an age spread (i.e., second-generation stars) in 
young massive star clusters. In fact, one of the prediction of 
the “in situ” scenario is that eMSTO should be observed also 
in young star clusters if they have the adequate properties 
in terms of mass and escape velocity. Young star clusters 
must be massive enough (i.e., M > 10® Mq) in order to 
have deep enough potentials and high enough escape veloc¬ 
ities to retain material lost due to stellar evolution of the 
first generation stars and to be able to accrete new material 
from the pristine gas still present in the surroundings. Un¬ 


fortunately, these strict constraints render young star clus¬ 
ters in which the presence of the eMSTO phenomenon is 
expected and verifiable very scarce. Glia identified the rel¬ 
atively young (~ 300 Myr) massive star cluster NGC 1856 in 
the LMC as a promising c a ndida te. Following this sugges¬ 
tion, IBastian fc Silva-Villal (l2013l j analyzed archival Wide 
Field & Planetary Camera 2 (WFPC2) images of this clus¬ 
ter, finding no evidence for age spreads larger than ~ 35 
Myr, which was interpreted by the authors as a suggestion 
that the eMSTO feature in intermediate-age clusters can- 
not be due to age spreads. However, the data analyzed by 
IBastian fc Silva-Villal (120131 ') seem to have reached only ~ 2 
mag beyond the MSTO, whose shape is vertical for the filters 
used in their study (F450W and F555W). This is not deep 
enough to sample a significant portion of the MS of the clus¬ 
ter, especially its curvature at lower stellar masses, which is 
important to constrain several important parameters of the 
cluster such as age, metallicity, distance, and reddening. In 
principle, this could have prevented the identification of fea¬ 
tures that might point to the presence of second-generation 
stars. 

With this in mind, we present an analysis of new HST 
WFC3 photometry of NGC 1856. These data allowed us to 
sample the cluster population down to ~ 8 mag beyond the 
MSTO, obtaining deep CMDs that permit a thorough ex¬ 
amination of the MSTO morphology and the nature of the 
cluster. We compare the observed CMD with Monte Carlo 
simulations in order to quantify whether it can be repro¬ 
duced by one simple stellar population (SSP) or whether it 
is best reproduced by a range of ages. We also study the evo¬ 
lution of the cluster mass and escape velocity from an age of 
10 Myr to its current age, to verify whether the cluster has 
the right properties to form and retain a second generation 
of stars. This analysis allows us to reveal new findings on 
the star formation history of the cluster NGC 1856 in the 
context of the multiple population scenario. 

The remainder of the paper is organized as follows: 
Section presents the observations and data reduction. In 
Section [3] we present the observed CMD, we describe the 
technique applied to correct the CMD due to the presence 
of differential reddening effects and we derive the best-fit 
isochrones. In Section[3]we compare the observed CMD with 
Monte Carlo simulations, both using one SSP and a combi¬ 
nation of two SSPs with different ages. We derive pseudo¬ 
age distributions for the observed and simulated CMDs and 
compare them. In Section we investigate the presence of 
ongoing star formation in the cluster, looking for PMS stars 
of the second generation. In Section[^we test the prediction 
of the stellar rotation scenario, while Section [7] presents the 
physical and dynamical properties of the cluster, deriving 
the evolution of the cluster escape velocity as a function of 
age. Finally, in Section [S] we present and discuss our main 
results. 


2 OBSERVATIONS AND DATA REDUCTION 

NGC 1856 was observed with HST on 2013 November 12 
using the UVIS channel of the WFC3 as part of the HST 
program 13011 (PI: T. H. Puzia). The cluster was centered 
on one of the two CCD chips of the WFC3 UVIS camera, 
so that the observations cover enough radial extent to study 
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Figure 1. Completeness fraction as function of magnitude 
F438W and distance from the cluster center. The open circles 
and so lid lines rep resent data within the core radius based 
on the iKind lll962r) model fit derived as described in Section 17.11 
open squares and dashed line represent data in the range between 
Ta and 2 xrc and the open triangles and dotted lines represent 
data outside 2 xrc. 

variations within cluster radius and to avoid the loss of the 
central region of the cluster due to the CCD chip gap. The 
cluster was observed in four filters, namely F438W, F555W, 
F658N, and F814W. Two long exposures were taken in each 
of the four filters: their exposure times were 430 s (F438W), 
350 s (F555W), 450 s (F814W) and 940 s (F658N). In ad¬ 
dition, we took three short exposures in the F438W, F656N 
and F814W filters (185 s, 735 s and 51 s, respectively), to 
avoid saturation of the brightest RGB and AGB stars. The 
two long exposures in each filter were spatially offset from 
each other by 2'.'401 in a direction -|-85?76 with respect to 
the positive X-axis of the detector, in order to move across 
the gap between the two CGD chips, as well to simplify 
the identification and removal of hot pixels and cosmic rays. 
In addition to the WFC3/UVIS observations, we used the 
Wide Field Camera (WFC) of ACS in parallel to obtain im¬ 
ages « 6' from the cluster center. These ACS images have 
been taken with the F435W, F656N and F814W filters and 
provide valuable information of the stellar content and star 
formation history in the underlying LMC field, permitting 
us to establish in detail the field star contamination fraction 
in each region of the CMDs. 

To reduce the ima ges we followed the method described 
in iKalirai et al.l l|2012l l. Briefly, we started from the jit files 
provided by the HST pipeline, which constitute the bias- 
corrected, dark-subtracted and flat-fielded images. After 
correcting the jit files for charge transfer inefficiency us¬ 
ing the dedicated CTE correction softwa rj^, we generated 
distortion-free images using MultiDrizzle dFruchter fc Hoo^ 

^ http://www.stsci.edu/hst/wfc3/tools/cte_tools 


R 

(1) 

F438W 

(2) 

F555W 

(3) 

F658N 

(4) 

F814W 

(5) 

R< Rc 

22.68 

22.50 

21.56 

22.06 

Rc < R <2- Rc 

24.78 

24.52 

22.33 

23.88 

R> 2-Rc 

26.55 

26.48 

22.60 

25.76 


Table 1.50% completeness limit for each bands in three different 
intervals in radius. (1): Radius interval (2-5): magnitudes in each 
band corresponding to 50% completeness fraction. 


Il997h and we calculated the transformation between the in¬ 
dividually drizzled images in each filter, linking them to 
a reference frame (i.e., the first exposure). Through these 
transformations we obtain an alignment of the individual 
images to better than 0.02 pixels. After flagging and reject¬ 
ing bad pixels and cosmic rays from the input images, we 
created a final image for each filter, combining the input 
undistorted and aligned frames. The final stacked images 
were generated at the native resolution of the WFC3/UVIS 
and ACS/WFG (i.e., 0'.'040 pixel”^ and 0'.'049 pixel”re¬ 
spectively). 

To perform the stellar photometry, we used the stand¬ 
alone versions of the DAOPHOT-II and ALLSTAR point 
spread function (PSF) fitting programs llStetsorJI 19871 . [l993 ~l 
on the stacked images. To obtain the final catalog we first 
performed the aperture photometry on all the sources that 
are at least 3a above the local sky, then we derived a PSF 
from ~ 1000 bright isolated stars in the field, and finally we 
applied the retrieved PSF to all the sources detected through 
the aperture photometry. We retained in the final catalogs 
only the sources that were iteratively matched between the 
two images and we cleaned them eliminating background 
galaxies and spurious detections by means of and sharp¬ 
ness cuts from the PSF fitting. 

Photometric calibration has been performed using a 
sample of bright isolated stars to transform the instru¬ 
mental PSF-fitted magnitudes to a fixed aperture of 10 
pixels (0'.'40 for WFC3/UVIS, 0'.'49 for ACS/WFG). We 
then transformed the magnitudes into the VEGAMAG sys¬ 
tem by adopting the relevant synthetic zero points for the 
WEC3/UVIS and ACS/WFG filters. Positions and magni¬ 
tudes, with the associated errors, for the first ten objects in 
our final catalog are reported in Table|XT]in the Appendix m 

We characterized the completeness as well as the pho¬ 
tometric error distribution of the final photometry by per¬ 
forming artificial star tests, using the standard technique 
of adding artificial stars to the images and running them 
through the photometric routines that were applied to the 
drizzled images using identical criteria. We added to each 
image a total of nearly 700000 artificial stars. To not induce 
incompleteness due to crowding in the tests themselves, the 
fraction of stars injected into a given individual image was 
set to be ~ 5% of the total number of stars in the final cat¬ 
alogs. The overall distribution of the inserted artificial stars 
was chosen to follow that of the stars in the image. We dis¬ 
tributed the artificial stars in magnitude in order to repro¬ 
duce a luminosity function similar to the observed one and 
with a colour distribution that span the full colour ranges 
found in the CMDs. After inserting the artificial stars, we 
applied to each image the photometry procedures described 


























































4 M. Correnti et al. 



Figure 2. Observed F438W vs. F438W — F814W CMD for all 
the stars inside the core radius, based on the iKin idi^) model 
fit derived as described in Section 17.11 Contamination from the 
underlying LMC field population has been derived from a region 
near the corner of the image, with the same surface area adopted 
for the cluster stars, and superposed on the cluster CMD (green 
open squares). Magnitude and colour errors, derived using the 
photometric distribution of our artificial stars within the core 
radius are shown in the left side of the CMD. The reddening 
vector is also shown for Ay = 0.4. 

above. The stars were recovered blindly and automatically 
cross-matched to the input starlist containing actual posi¬ 
tions and fluxes; they were considered recovered if the input 
and output magnitudes agreed to within 0.75 mag in both 
filters. Finally, we assigned a completeness fraction to each 
individual star in a given CMD as a function of its magni¬ 
tude and distance from the cluster center. The completeness 
fraction of stars as a function of the magnitude F438W and 
distance from the cluster center, in three different intervals 
in radius, are shown in Figure [T] The magnitudes corre¬ 
sponding to the 50% completeness fraction in each band, 
for the same intervals, are reported in Table [T] 


3 COLOUR-MAGNITUDE DIAGRAM 
ANALYSIS 

3.1 A wide main sequence turn-off region 

Figure [3 shows the F438W - F814W vs. F438W CMD of 
NGC 1856, plotting only the stars within the core radius, 
based on the iKinj 1 19621 ) model fit derived as described in 
Section rm The observed CMD presents some interesting 
features. While the lower part of the MS (i.e., below the 
“kink” at F438W ~ 21.2) is quite narrow and well defined, 
the MSTO region (with F438W < 19.5), appears relatively 
wide when compared with the photometric errors. 

To determine whether the observed broadening of the 
MSTO region is a “real” feature in the CMD, we first test 
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Figure 3. CMDs for the 9 subregions in which we divided the 
cluster field, derive d as described in Sec tion [3] Best-fit isochrones 
(black lines) from iMarigo et al.l feOPSl l are superposed in each 
CMD, along with the derived visual extinction Ay. The top leg¬ 
end shows the age, distance modulus (m — M)o and metallicity 
adopted for all the isochrones. 

whether it could be caused by the following effects: contam¬ 
ination by the underlying field population, poor photometry 
(i.e.; large photometric uncertainties), the presence of differ¬ 
ential reddening, and/or a significant binary fraction. 

To assess the level of contamination of the underlying 
LMC field population, we selected a region near the corner 
of the image with the same surface area adopted for the 
cluster stars. Stars located in this background region have 
been superposed on the cluster’s CMD (shown as green open 
squares in Figure [2]) • The contamination is mainly confined 
to the lower (faint) part of the MS and does not affect the 
MSTO and the Red Clump (RC) in any significant fashion; 
thus we can conclude that underlying LMC field stars do 
not cause the observed shape of the MSTO or RC regions. 

Photometric uncertainties have been derived from the 
artificial stars test. Magnitude and color errors are shown 
in Figure [3 photometric errors at the MSTO level are be¬ 
tween 0.02 and 0.04 mag, far too small to account for the 
broadening of the MSTO. 

3.2 The impacts of differential reddening and 
binary fractions 

To check for the presence of differential reddening in the 
cluster field we adopted the following approach. First of all, 
we selected a square region that circumscribes the area de¬ 
fined by the core radius (i.e.; the length of the side of the 
square is equal to twice the core radius), we divided this re¬ 
gion into 9 equal-sized squares, and we compared the CMDs 
derived in each of them. The number of equal-sized region 
has been arbitrarily choosen in order to have a good sam¬ 
pling of the different region mantaining a statistically sig- 
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Figure 4. Spatial distribution of the differential reddening in 
NGC 1856. The color of each star represents the final AE{B — I) 
derived from the differential reddening correction inside the core 
radius region. The color coding is shown at the top. 


nificant number of objects in each field. This approach has 
been adopted only to check if differential reddening effects 
are present in the cluster field and to have a rongh estimate 
of the reddening variation. Figure [3] shows the CMDs for 
the 9 subregions; to yield a preliminary hint on the amount 
of differential r eddening, we superim posed on each CMD an 
isochrone from iMarigo et ahl (120081 ') for which age, distance 
modulus and metallicity are fixed (see top legend of Fig¬ 
ure 0. while the visual extinction Ay is left free to vary. 
Figure |3] clearly shows that the best fit is achieved in each 
CMD using a different value of the visual extinction Ay 
(derived values of Ay are reported in each field), confirming 
that a differential reddening is present in the cluster field, 
with variations of the order of ~ ± 0.15 mag with respect 
to the mean reddening value. 

Once verified that the differential reddening is truly 
present in the cluster field, we corrected each star in the 
CMD using t he following approa ch (for a more detailed de¬ 
scription, see lMilone et al.ll2012l '). Briefly, we first defined a 
photometric reference frame in which the X axis is parallel 
to the reddening line. In this reference system, it is much 
easier to determine reddening difference rather than in the 
original colour-magnitude plane, where the reddening vector 
is an oblique line. To do this, we arbitrarily defined an origin 
O, then translated the CMD such that the origin of the new 
reference frame corresponds to O and then we rotated the 
CMD counterclockwise by an angle defined by the equation: 


6 = arctan 


(- 


Apissw 




VMf438W — Af&UW J 


( 1 ) 


where Afassw and Afsi 4 w are the appropr iate extinction 
coefficients for the UVIS WFC3 filters. Using ICardelli et al.l 
lll989ll and lO’Donn^ (Il994h extinction curve and adopt¬ 


ing Ry = 3.1, their values are ^f 438W = 1.331'Av and 
^F8i4W = 0.610'^v, respectively. At this point, we gen¬ 
erated a fiducial line using only MS stars. We divided the 
sample into bins of fixed “magnitude” and we calculated the 
medians in the X- and Y-axis. The use of the median allows 
us to minimize the influence of outliers such as binary stars 
or stars with poor photometry left in the sample after the 
selection. Among the MS stars, we selected a subsample lo¬ 
cated in the region where the reddening line defines a wide 
angle with the fiducial line, so that the shift in colour can be 
safely interpreted as an effect of the differential reddening. 
For this reason, we limited our sample to the central portion 
of the cluster MS, excluding the upper part near the MSTO 
and the lower portion, fainter than the magnitude at which 
the MS starts to bend in a direction parallel to the redden¬ 
ing line. For each selected star, we calculated the distance 
from the fiducial line along the reddening direction (i.e. A 
X) and we used these stars as reference stars to estimate the 
reddening variations associated to each star in the CMD. Fi¬ 
nally, to correct them, we selected the nearest 30 reference 
stars and we calculate the median A X that is assumed as 
the reddening correction for that particular star (to derive 
the differential reddening suffered by the reference stars, we 
excluded that star in the computation of the median A X). 
We note that the median correlation length introduced with 
this approach is of the order of ~ 20 px (corresponding to 
0.2 pc). After the median A X values have been subtracted 
from the X-axis value of each star in the rotated CMD, we 
obtained an improved diagram that we used to derive a more 
accurate selection of MS reference stars and a more precise 
fiducial line. We then applied again the described procedure, 
iterating the process until it converges (in this case, a couple 
of iterations were sufficient). At this point, the corrected X 
and Y “magnitude” were converted back to the F438W and 
F814W magnitudes. 

Figure|4]shows the spatial distribution of the differential 
reddening inside the core radius of NGC 1856; each star 
is reported with a different color, depending on the final 
Ai?(F438W —F814W) applied to correct it, derived with the 
method described above. In this context, we observe that the 
CMD in the central panel of Figure 0 which represents the 
innermost region of the cluster, shows a quite wide MSTO, 
compared to the other panels. 

Conversely, in the corresponding spatial region from 
which the CMD is drawn (i.e., the central part in the map 
of Figure 0 the derived reddening is quite uniform, with 
Ai?(F438W — F814W) ~ 0.04 (corresponding to AAy ~ 
0.055). Moreover, it is worth to note that the reddening cor¬ 
rection reaches its highest accuracy in this inner region, due 
to the high surface number density of stars (i.e., smaller dis¬ 
tance of the star to be corrected from the reference stars 
used to derive AAy, see also Section EJ. 

We acknowledge that our differential reddening correc¬ 
tion might be somewhat ineffective for a small fraction of 
stars, depending on their location in the field, e.g., in the 
outer regions, where the distance between stars is generally 
larger than in the inner regions. However, the number of 
stars in such regions is by nature relatively low (see, e.g., 
the CMDs in the corner panels of Figure 0. Thus, only a 
small fraction of MSTO stars can suffer from a less accu¬ 
rate correction, leaving the global morphology of the MSTO 
region unaltered. 
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F438W - F814W 


Figure 5. Differential-reddening-corrected F438W - F814W versus F438W CMD for NGC 1856. Best-fit isochrones from iMarigo et aTI 
ll2()08l i. for the minimum (red solid line) and maximum (blue dashed line) ages (300 and 410 Myr, respectively) that can be accounted 
for by the data are superposed on the cluster CMD, along with the derived distance modulus (m — M)o and visual extinction Ay. The 
parallelogram box used to select MSTO stars in order to derive the pseudo-age distribution, as described in Section |4] is also shown. 


The differential-reddening-corrected (DRC) CMD is 
shown in Figure [5] The correction caused the lower part 
of the MS to be even more defined and narrow than it 
was in the observed CMD. Conversely, the MSTO region 
and upper MS is still fairly wide, although it has narrowed 
down slightly. This is illustrated in more detail in Figure [S] 
which shows the distribution of F438W — F814W colours 
in the magnitude range 18.7 ^ F438W ^ 19.mag for 
the observed CMD and that obtained for the DRC CMD 
(red dashed and black solid lines, respectively). Note that 
the global shape and width of the two distributions is very 
similar, except for the fact that the observed one has more 
extended outer wings than the one derived from the DRC 
CMD. These extended wings are caused by stars that are 
located in the region where the differential reddening effects 
are stronger and then have been shifted from their original 
position in the CMD. The fact that the shape of the two 
distributions is very similar within their FWHM values sug¬ 
gests that the observed broadening in the MSTO region is an 
intrinsic feature, rather than due to differential reddening. 

Given the above, we judge it very unlikely that differ- 


^ This magnitude range was selected in this context because it 
reflects the part of the upper MS for which the binary sequence 
is merged in with the single star sequence (see Fig. (7)1. Hence, our 
result is robust against any given binary fraction. 


ential reddening effects or binary stars can explain the ob¬ 
served broadening of the MSTO region. 


3.3 Isochrone fitting 

Isochrone fitting was done following the methods described 
in detail in iGoudfrooii et al.l (l2011bl . hereafter Glib) for the 
case of star clusters with ages that are too young to have a 
well-developed red giant branch. Briefly, age and metallicity 
were derived using the observed differences in magnitude 
and colour between the MSTO and the RG; we selected all 
the isochrones for which the values of those parameters lie 
within 2(7 of the uncertainties of those parameters as derived 
from the GMDs. For the set of isochrones that satisfied our 
selection criteria (5-10 isochrones), we found the best-fit 
values for distance modulus and reddening by means of a 
least-square fitting program to the magnitudes and colors of 
the MSTO and RG. Finally, we overplotted the isochrones 
onto the CMDs and selected the best-fitting ones by means 
of a visual examination. 

In this cont ext, we superposed o n the cluster GMD two 
isochrones from iMarigo et al.l ll2008l) of different ages (300 
Myr and 410 Myr, respectively), shown in Figure[5]with dif¬ 
ferent colours and line styles (red solid line for the younger 
isochrone and blue dashed line for the older one, respec¬ 
tively). The isochrones ages were chosen to match the min- 
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Figure 6. Comparison between the colour distributions obtained 
from the differential-reddening-corrected CMD (black solid line) 
and the observed CMD (red dashed line). 


imum and maximum age that can be accounted for by the 
data and isochrone fitting performed as described above. 

Taking the results shown in Sections 13.11 - 13.31 at face 
value, it seems that the morphology of the MSTO is not 
very well reproduced by a SSP, and that a spread in age 
of the order of ~ 100 Myr may constitute a better fit to 
the data. This is further addressed below (Section 14.11 and 
Section IT^ . 


4 MONTE CARLO SIMULATIONS 

To further determine whether or not a single population can 
reproduce the observed CMD, we conducted Monte Carlo 
simulations of a synthetic cluster with the properties im¬ 
plied by the isochrone fitting (see Gila and G14 for a de¬ 
tailed description of the method adopted to produce these 
simulations). Briefly, to simulate a SSP with a given age 
and chemical composition, we populated an isochrone with 
stars randomly drawn using a Salpeter mass function and 
normalized to the observed (completeness-corrected) num¬ 
ber of stars. To a fraction of these sample stars, we added 
a component of unresolved binary stars derived from the 
same mass function, using a flat distribution of primary-to- 
secondary mass ratios. The binary fraction was estimated 
by using as a template the part of the observed MS between 
the MSTO region and the TO of the background population 
in the magnitude range 22.5 < F438VK < 19.5 (see Fig. [2]). 
We derived that the value of the binary fraction that best 
fit the data is of ~ 25%. We estimated that the internal 
systematic uncertainty in the binary fraction is ~ 5%; for 
the purposes of this work, the results do not change signifi¬ 
cantly within 10% of the binary fraction. Finally, we added 
photometric errors that were derived using our artificial star 
tests. We derived three different sets of synthetic CMDs: one 


with a single SSP of age 300 Myr, and two combining two 
SSPs of different ages (one with SSPs with ages of 300 Myr 
and 380 Myr, and the other with ages of 300 Myr and 410 
Myr, respectively). In the last two cases, we derived a set of 
simulations in which we changed the ratio of the number of 
stars in the two populations, from 10% to 90% of stars in 
the younger SSP, with a 5% step increase in each run. 

To compare in detail the observed MSTO region with 
the simulated ones, we created “pseudo-age” distributions 
(see Gila for a detailed description). Briefly, pseudo-age 
distributions are derived by constructing a parallelogram 
across the MSTO, with one axis approximately parallel 
to the isochrones and the other approximately parallel to 
them (the adopted parallelogram is shown in Figure[5]). The 
(F438W — F814W, F438W) coordinates of the stars within 
the parallelogram are then transformed into the reference co¬ 
ordinate frame defined by the two axis of the parallelogram 
and the same procedure is applied to the isochrones tables to 
set the age scale along this vector. The pseudo-age distribu¬ 
tions are calculated usin g the non-parame tric Epanechnikov- 
kernel density function ( Silvermcinlll986l ). in order to avoid 
possible biases that can arise if fixed bin widths are used. 
Finally, a polynomial least-squares fit between age and the 
coordinate in the direction perpendicular to the isochrones 
yields the final pseudo-age distributions. The described pro¬ 
cedure is applied both to the observed and the simulated 
CMD. In the following sections, we describe the results ob¬ 
tained from the comparison of the observed and simulated 
pseudo-age distributions, for the synthetic CMDs obtained 
with a single SSP and for those obtained from the combina¬ 
tion of two SSPs of different ages. 


4.1 Comparison with the synthetic CMD of a SSP 

The top panel of Figure [7] shows the comparison between 
the observed and simulated CMD, the latter obtained from 
a single SSP with an age of 300 Myr. Overall, the SSP simu¬ 
lation reproduces the CMD features quite well in the fainter 
portion of the cluster MS (i.e.; for F438W > 20.5 mag). Con¬ 
versely, we note that the MSTO regior0of the simulated SSP 
does not reproduce the observed MSTO region very well, in 
that the latter is wider than the former. Furthermore, the 
faint end of the RC also does not seem to be reproduced well 
by the simulated SSP in that the simulated luminosity func¬ 
tion of RC stars peaks more strongly at bright magnitudes 
than the observed one. Note that this is consistent with a 
fraction of stars in NGC 1856 having “older” ages than that 
of the best-fit SSP (see for reference the shape and location 
of the RC for the older isochrone in Figure [5]). This differ¬ 
ence is even more evident in a comparison of the observed 
and simulated pseudo-age distributions (see bottom panel 
of Figure 0. Indeed, the observed pseudo-age distribution 
(shown as a black solid line in Figure0 reaches significantly 
older ages than the simulated one (red dashed line in Fig¬ 
ure 0, while the two distributions are very similar to each 
other in the left (“young”) half of the respective profiles. 


® i.e., the part of the MS with F438W < 19.2 mag, where 
the binary sequence has joined the single star sequence in 
F438W - F814W 
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Figure 7. Top panel: comparison between DRC (gray dots) and 
simulated CMDs (red circles, single stars, magenta square, binary 
stars, respectively). The simulated CMD is obtained from Monte 
Carlo simulation of a SSP with an age of 300 Myr. The paral¬ 
lelogram box used to select MSTO stars and the fraction of the 
adopted binary stars are also reported. Bottom panel: Pseudo-age 
distribution for the MSTO region of the DRC (black solid line) 
and simulated (red dashed line) CMDs. 


To quantify the difference between the pseudo-age dis¬ 
tribution of the cluster data and that of the SSP simulation 
in term of intrinsic MSTO width of the cluster, we measured 
the widths of the two sets of distributions at 50% of their 
maximum values (hereafter called FWHM), using quadratic 
interpolation. The intrinsic pseudo-age range of the cluster 
is then estimated by subtracting the simulation width in 
quadrature: 


Figure 8. Comparison between the observed pseudo-age distribu¬ 
tion (from the DRC CMD, black solid line) and the ones obtained 
from the combination of two SSPs with different ages (300 Myr 
and 410 Myr, top panels; 300 Myr and 380 Myr, bottom panels) 
shown as red dashed lines. In each panel we report the adopted 
ages of the two SSPs and the number ratio between the younger 
SSP and the old one. Two-tail p values obtained for the K-S test 
are also shown. 


FWHMmsto = {FWHMI^, - FWHMlsvf^^ 


( 2 ) 


where the “obs” subscript indicates a measurement on 
the DRC CMD and the “SSP” subscript indicates a 
measurement on the simulated CMD for the SSP. With 
FWFIMohs = 269 Myr and FWFIMssp = 198 Myr, we ob¬ 
tain FWHM MSTO ~ 182 Myr, which is similar to the value 
of FWHMssp- This suggests that NGC 1856 may host two 
SSPs separated in age by about FWHM msto /‘2 ~90 Myr 
(see Sect. 14.21 belowl. 

In this context, we adopted equation [2] also to derive 
the value of reddening variation necessary to reproduce the 
observed pseudo-age distributions assuming that a cluster 
is formed by a SSP. Instead of age, we derived the FWHM 
of the observed and simulated pseudo-age distributions as a 
function of F438W - F814W; using the appropriate extinc¬ 
tion coefficient for the UVIS WFC3 filters we transformed 
the “A color” in AAy obtaining the following value: AAy 
= 0.27. The derived AAv is significantly larger than the 
reddening variations observed in the spatial differential red¬ 
dening map in the innermost region of the cluster. More¬ 
over, we note that due to the uniform shape of the observed 
pseudo-age distribution, this large reddening should affect 
a significant number of stars, making the possibility that 
the reddening can mimic the presence of a multiple popu¬ 
lation even less plausible. Therefore, taking these results at 
face value, they seem to suggest that a SSP is not able to 
reproduce the observed morphology of the MSTO region. 
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Figure 9. Left panels: DRC CMDs zoomed in the upper po rtion of the MS with superimposed isochrones (red solid and blue dashed 
lines for the young and the old isochrones, respectively) from iMarigo et alj 120081 ) along with the adopted ages (left panel: 300 and 410 
Myr. Right panel: 300 and 380 Myr). Dotted lines represent the magnitude and colour cuts adopted to select RC stars, whereas the 
dotdashed lines represent the magnitude cut, used to divide the two parts of the RC. Right panels: ratios between the stars in the upper 
portion of the RC and their total number as a function of the ratio of the young SSP (top panel: SSPs with ages of 300 and 410 Myr; 
bottom panel: SSPs with ages of 300 Myr and 410 Myr). The ratio obtained in the observed DRC CMD is reported as a dashed line in 
both panels. The typical uncertainty (^ 10%) in the calculated number ratio is also shown. 


4.2 Comparison with a synthetic CMD of two 
SSPs of different age 

As stated in Section IH we also simulated synthetic CMDs 
combining two SSPs with different ages. One set is obtained 
using SSPs with ages of 300 Myr and 410 Myr and the other 
with ages of 300 Myr and 380 Myr. For both cases, we de¬ 
rived a set of simulations in which we varied the ratio of 
the number of stars in the young and the old SSP, from 
10% to 90% in each of them, with a step of 5%. For each 
simulated CMD we obtained the corresponding pseudo-age 
distribution that we compared with the observed one. 

The fact that the stars from the two SSPs are mixed 
together in the parallelogram used to derive the pseudo¬ 
age distribution, combined with our purely statistical ap¬ 
proach, cause an increase in the “free parameters” and in¬ 
troduces a sort of degeneracy in the derived results. This is 
reflected in the fact that we obtained more than one sim¬ 
ulated pseudo-age distribution that reasonably reproduces 
the observed one. These “best-fitting” pseudo-age distribu¬ 
tions are shown in Figure [S] In detail, the top panels show 
the best fit achieved from the SSPs with ages of 300 Myr 
and 410 Myr (with mass fractions of “young” stars of 40% 
and 50%, respectively), whereas in the bottom panels we 
show the ones obtained from the SSPs with ages of 300 Myr 
and 380 Myr (with mass fractions of young stars of 40% and 
55%, respectively) The ages of the two SSPs and the ratio 
between them have been reported in each panel. We note 
that the best-fit mass ratios of the young to the old popu¬ 
lation seem to be around 1:1, consistent with the observed 
pseudo-age distribution, where we see a hint of two peaks 


with a similar maximum value and hence a likely number 
ratio of ~ 1:1 for young vs. old stars in the cluster. 

To constrain which combination of two SSPs pro¬ 
vides the best solution to reproduce the observed CMD of 
NGC 1856, we performed two-sample Kolmogorov-Smirnov 
(K-S) tests. The two-tailed p values are mentioned in Fig¬ 
ure [S] for each pseudo-age distribution. Taking these results 
at face value, it seems that the combination of SSPs with 
ages of 300 and 380 Myr provide a better solution with re¬ 
spect to the SSPs with ages of 300 and 410 Myr. To further 
constrain the two-SSP fits, we pointed our attention to the 
distribution of RC stars, adopting the following approach. 
The left panels of Figure [^show the DRC CMDs zoomed in 
to the RC and the up per portion of t he M S, with superim¬ 
posed isochrones from iMarigo et ^ (l2008l ). with the same 
ages as those adopted in the SSPs (left sub-panel: 300 and 
410 Myr; right sub-panel: 300 and 380 Myr, respectively). 
First, we selected cluster RC stars by applying colour and 
magnitude cuts (shown as dotted lines in the left panels of 
Figure^. Then, we divide RC stars in two parts, applying a 
magnitude threshold (shown as dotdashed lines in Figure[9ll 
coinciding with the brightest point of the RC of the “old” 
isochrone (i.e., F438W = 19.3 mag for the 410 Myr isochrone 
and F438W = 19.15 mag for the 380 Myr isochrone). Finally, 
we counted the stars in the upper and lower portion of the 
RC, both in the observed data and in each synthetic CMD 
obtained from the associated Monte Carlo simulations. 

The ratios between the stars in the upper portion of the 
RC and their total number as a function of the mass ratio 
of the young SSP in each simulation are shown in the right 
panels of Figure [9] (for SSPs with ages of 300 Myr and 410 
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F438W - FdUW F438W - F814W 


Figure 10. Top panels: F555W - F814W vs F555W - F658N 
colour-colour diagram for the stars inside the core radius (left 
panel) and in the control field (right panel), representing the con¬ 
tamination from the background LMC stars (same field adopted 
in Figure [a. The dashed line represents the median F555W - 
F658N colour, representative of stars with no Ha excess. Stars 
with a F555W - F658N colour at least bcr above that reference 
line, where a is the uncertainty on the F555W - F658N colour 
of the star, are reported as open red squares. Those among them 
with EWho. > lOA, where the equivalent widt h of the Hg emis¬ 
sion li ne {EWho. ) calculated using the Eq. 4 of lPe Marchi et al.l 
ll2010l'l are reported as solid blue circles. Bottom Panels: CMDs 
of the cluster and of the control field with overplotted stars se¬ 
lected from the F555W - F814W vs F555W - F658N colour-colour 
diagram (solid blue circles). 

Myr in the top panel and for SSPs with ages of 300 Myr and 
380 Myr in the bottom panel, respectively). The ratios from 
the DRC CMD, obtained with the different magnitude cuts, 
are shown as dashed lines. The typical uncertainty on the 
mass ratios has been estimated to be of order 10 %, shown 
as an error bar in both panels. Taking these results at face 
value, it seems that the best match is obtained from the 
synthetic CMD obtained from the SSPs with ages of 300 
Myr and 380 Myr and with a fraction of young stars of 
0.55 ±0.10 (cf. bottom right panel in Figure [ 8 ll. 

Globally, the results presented in Section 14.11 and EJ 
show that a single SSP fails in reproducing the observed 
pseudo-age distribution while a combination of two SSPs 
with different ages can provide a very good fit. This seems 
to indicate that a spread in age may be present in the 
cluster and can not be cate gorically excluded as argued by 
iBastian fc Silva-Villal (l2013l ). 

5 CONSTRAINTS ON ONGOING STAR 

FORMATION ACTIVITY 

To study whether NGC 1856 presently hosts ongoing star 
formation, we performed a search for PMS stars by means of 


narrow-band imaging of the Hq line. Ha emission is a good 
indicator of the PMS stage: the presence of strong Ha emis¬ 
sion {EWho > 10 A) in young stellar objects is normally 
interpreted as a signature of the mass accretion process onto 
the surface of the object that requires the presence of an 
inner disk (see lFeigelson fc Montmerlgl 19991 : 1 White fc Basril 
l2003l . and re ference therein). To do t his, we used the method 
described in lPe Marchi et ^ ll2010l) . which combines V and 
I broad-band photometry with narrow-band Ha imaging to 
identify all the stars with excess Ha emission and to mea¬ 
sure their Ha luminosity. 

In the top panels of Figure [in] we show the 
F555W — F814W vs. F555W — F658N colour-colour dia¬ 
gram for the stars inside the cluster core radius (left top 
panel) and for the stars in the “control field”, the same 
field we used to derive the contamination by the back¬ 
ground LMC population in Figure [2] (right top panel). 
We used the median F555W —F658N colour of stars with 
small (< 0.05 mag) photometric uncertainties in each of the 
three F555W, F814W, and F658N bands, as a function of 
F555W — F814W, to define the reference template with re¬ 
spect to which excess Ha emission is identified (shown by 
the dashed line in Figure nni- We selected a first candi¬ 
date sample of stars with Ha excess emission by consid¬ 
ering all those with a F555W — F658N color at least 5a 
above that reference line, where a is the uncertainty on 
the F555W — F658N colonr of the star in question (shown 
as open red squares in Figure fTOl) . Then we calculated the 
equivalent width of the Ha emission line {EWhc) from the 
measured color e xcess , using the following equation from 
IPe Marchi et ah lim^: 

Weg{Ha) = RW X [1 - ( 3 ) 

where RW is the rectangular width of the filter (similar in 
definition to the equivalent width of the line), which depends 
on the characteristic of the filter (for the adopted F658 N fil¬ 
ter, Rw = 17.68; see Table 4 in iDe Marchi etliDbOld ). Fi¬ 
nally, we considered as bona-fi de candidate PMS sta rs those 
objects with EWhc > 10 A ll White fc Basril [20031 ) . shown 
as solid blue circles in Figure HOI The number of objects 
that show excess Ha emission is very low and is roughly 
the same in both fields. In the bottom panels of Figure fTCTl 
we overplotted these objects (shown as solid blue circles) on 
the F438W - F814W vs. F438W GMDs of the cluster and 
the control field; the majority of them are located on the 
fainter part of the GMD, where the photometric errors are 
larger and the contamination by background LMC stars is 
higher. In fact the number of objects that show an Ha ex¬ 
cess is comparable within the errors in the two fields below 
the TO of the background stellar population (i.e. F438W 
< 23 mag). Hence, the number and the location of these 
objects in the CMDs suggest that they can be considered 
spurious detections. For what concerns the handful of ob¬ 
jects at brighter magnitudes (i.e. F438W > 22 mag), where 
the photometric errors are smaller and the contamination 
is lower, we hypothesize that these can be stars in a binary 
system in which mass transfer between the primary and the 
secondary star is occurring. We exclude the possibility that 
these objects are true PMS stars due to the fact that, in this 
case, we would have observed a significant number of these 
objects at fainter magnitudes (i.e. lower masses and hence 
longer PMS lifetime). 
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Figure 11. Co mparison between th e non-rotating isochrone (blue 
solid line), from iMarigo et ^ ||2008^ , and its rotating counterpart 
(rec dashed line), derived as described in Section|6l Isochrones are 
superposed on the cluster DRC CMD, along with adopted age, 
distance modulus (m — M)o and visual extinction Ay. 

In this context, it is worth to note that our magnitude 
detection limit for stars showing Ha excess is around F438W 
~ 26.5 mag, which, at the cluster age, corresponds to stars 
with masses of ~ 0.8 M©. The PMS lifetime for such stars is 
of the order of ~ 120 Myr, indicating that if star formation 
is occuring within the cluster, with our method we should 
be able to observe a significant number of PMS objects. 

Hence, these results indicate that the level of ongoing 
star formation activity is negligible in NGC 1856. In fact, 
the lack of low-mass stars showing an Ha excess indicates 
that the star formation in the cluster must have stopped 
at least ~ 120 Myr ago, in agreement with the conclusions 
derived from the pseudo-age distribution analysis. 


6 CONSTRAINTS TO THE STELLAR 

ROTATION SCENARIO IN PRODUCING 
EMSTO’S 

In the literature, some alternative explanations that do 
not invoke the presence of an extended star formation, 
have been proposed to explain the eMSTO phenomenon 
in intermediate-age star cluster s. In particular, the so 
called “stel l ar rot ation” scenario dBastian fc de Minkll200g : 
lYang et al.l l2013l ) suggest that eMSTOs can be explained 
by a spread in rotation velocity among turn-off stars. 
lYang et al.l ll2013l) calculated evolutionary tracks of non¬ 
rotating and rotating stars for three different initial stel¬ 
lar rotation periods (approximately 0.2, 0.3 and 0.4 times 
the Keplerian rotation of ZAMS stars), and for two differ¬ 
ent rotational mixing efficiencies (“normal”, fc = 0.03 and 
“enhanced”, fc = 0.20). From the isochrones, built from 
these tracks, they calculated the widths of the MSTO region 


caused by stellar rotation as a function of cluster age and 
translated t hem to age spreads. In particular, Figure 7 in 
lYang et al.l (l2013l ) shows the equivalent width of the MSTO 
of star clusters caused by rotation as a function of the clus¬ 
ter age, for different initial stellar rotation periods and ro¬ 
tational mixing efficiencies. At the age of NGC 1856, all the 
Yang et al. models with “normal” mixing efficiency (i.e.; 
fc = 0.03) indicate that rotation does not cause any ap¬ 
preciable spread. On the other hand, for their models with 
“enhanced” mixing efficiency (i.e.; fc = 0.20), the predic¬ 
tions vary from no spread for the model with the longest 
initial stellar rotation period to a maximum of ~ 100 Myr 
for the model with the shortest period (i.e., the model for the 
highest rotation velocity). In this case their model predicts 
a “negative” spread, in the sense that the rotating model 
exhibits a bluer MSTO with respect to their non-rotating 
counterpart, thus mimicking the presenc e of a younger po p- 
ulation (see top left panel of Figure 6 in lYang et al.ll201^ . 

To test the prediction of this rotating model in the spe- 
cific case of NGC 18 56, we used the 410 Myr isochrone from 
iMarigo et al.l ll2008l ) as a template for non-rotating stars 
(i.e., the same isochrone as that plotted in Figure[5ll, and we 
derived its rotating counterpart. To do so, we calculated the 
shifts in colour and magnitude betwee n the rotat i ng an d 
non-rotating isochrones in Figure 6 of lYang et al.l ll2013l) . 
we then transformed these s hifts in our photome tric system 
and we applied them to the lMarigo et all (l2008f) isochrone. 
Figure [TT] shows the DRC CMD of NGC 1856 with superim¬ 
posed non-rotating and rotating isochrones (blue solid and 
red dashed line, respectively). For the rotating isochrone, we 
plotted only the MS part, since the post-MS era of stellar 
evolution is not addressed by the lYang et al.l (l2013f) models. 
Overall, the rotating model plotted in Figure [TT] seems to 
reproduce the observed CMD reasonably well in the upper 
portion of the MS, whereas the bottom portion of the MS is 
not as well reproduced, the rotating model being somewhat 
too red with respect to the observed MS. The latter is due 
to the fact that the rotating model bends toward red colors 
(lower temperatures) for stellar masses Ai < 2.0 M©. 

In conclusion, the rotating model of lYang et al.l (l2013fl 
that involves “enhanced” rotational mixing efficiency seems 
to provide a satisfactory fit to the observed MS of NGC 1856, 
indicating that a spread in rotation velocity, beside an age 
spread, could be a possible cause of the observed broadening 
of the MSTO region. Ho wever, we note that the predictions 
of the lYang et al.l ll2013l) models are not consistent with the 
observations for intermediate-age (~ 1- 2 Gyr) star clusters 
featuring eMSTOs. This is especially the case for the rotating 
models involving enhanced rotational mixing efficiency (see 
discussion in G14, in particular their Figure 7). It should 
however also be recognized that the study of the creation of 
theoretical stellar tracks and isochrones for rotating stars at 
various stages of stellar evolution, rotation rates, and ages 
is still in relatively early stages. To date, no stellar rotation 
velocity measurements have yet been undertaken in young 
and intermediate-age star clusters in the Magellanic Clouds. 
We strongly encourage such studies in the future in order to 
provide fundamental improvements of our understanding of 
the possible relation between the eMSTO phenomenon and 
the effects of stellar rotation. 
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7 INSIGHTS FROM DYNAMICAL ANALYSIS 

As mentioned in the Introduction, a peculiar characteristic 
of eMSTOs in intermediate-age star clusters in the Magel¬ 
lanic Clouds is that they are not hosted by all such star 
clusters. To explain this phenomenon in the context of mul¬ 
tiple stellar populations, Gila proposed a scenario in which 
eMSTOs can only be hosted by clusters for which the escape 
velocity of the cluster is higher than the wind velocities of 
“polluter” stars thought to provide the material out of which 
the second stellar generation is formed, at the time such stars 
were present in the cluster (we refer to this as “the escape ve¬ 
locity threshold” scenario). This scenario was developed fur¬ 
ther by G14 who studied HST photometry of a sample of 18 
intermediate-age (1-2 Gyr) star clusters in the Magellanic 
Clouds that covered a variety of masses and core radii. They 
found that all the clusters showing an eMSTO feature escape 
velocities \4sc ^ 15 kms“^ out to ages of at least 100 Myr. 
This age is equivale nt to the lifetime of stars of ~ 5 Mq (e.g., 
iMarigo et al.ll2003) , and hence old enough for the slow winds 
of massive binary stars and IM-AGB stars of the first gener¬ 
ation to produce significant amounts of “polluted” material 
out of which second-generation stars may be formed. Fur¬ 
thermore, C14 showed that this threshold on 14sc is consis¬ 
tent with HST observations of four low-mass intermediate- 
age star clusters (~ 10^ Mq)\ the two clusters that host 
eMSTOs have 14sc ^ 15 kms“^ out to ages of ~ 100 Myr, 
whereas Ksc < 12 kms“^ at all ages for the two clusters 
not exhibiting eMSTOs. Hence, the critical escape velocity 
for a star cluster to be able to retain the material ejected 
by first-generation polluter stars seems to be in the range 
of 12- 15 kms“^. This threshold is consistent with wind ve¬ 
locities of IM-AGB stars and massive binary stars (see G14 
for a detailed discussion). Briefly, IM-A GB stars show wind 
veloc i ties in the range 1 2-18 km s~^ ilVassiliadis &: WoodI 
1 19931 : IZiilstra et al.l 1 19961 ). whereas observed velocities of 
ejecta of massive binary stars are in the range 15- 5(0 
kms~^ llSmith et al.l[200a. 120071') . 

With this in mind, following the results presented in 
the previous sections where the analysis of pseudo-age dis¬ 
tributions seems to suggest the presence of multiple stellar 
populations in NGC 1856, we determined the structural pa¬ 
rameters and the dynamical properties of the cluster from 
our new HST/WFC3 data. 



Figure 12. Radial surface number density profile of NGC 1856. 
Black points represent observed values. The dashed line represents 
the best-fit King model (cf. equation [4|l whose parameters are 
shown in the legend. Ellipticity and effective radius rg are also 
shown in the legend. The radius values have been converted to 
parsec by adopting the appropriate distance modulus. 









Figure 13. Cumulative completeness-corrected radial distribu¬ 
tion of bright versus faint stars. Red solid circles and solid line 
represent stars with F438W < 20.5 mag while black open circles 
and dotted line represent stars with 21.0 < F438W < 23.5 mag. 


7.1 Structural parameters 

We determined the radial surface number density dis- 
tribution of star s , foll owing the procedure described in 
iGoudfrooii et al.l ( 20091') . Briefly, we first determined the 
cluster centre to be at reference coordinate (xc, t/c) = (2986, 
1767) with an uncertainty of ± 5 pixels in either axis. In 
order to derive it, we first created a two-dimensional his¬ 
togram of the pixel coordinates adopting a bin size of 50 
X 50 pixels and then calculated the centre using a two- 
dimensional Gaussian fit to an image constructed from the 


^ We note that this measure have been derived from one sys¬ 
tem, RY Scuti, which is in a specific case of stable mass transfer, 
whereas most of the mass is ejected during the unstable phase. 
Lower velocities are expected during unstable mass transfer or 
the ejection of a common envelope. 


surface number density values in the aforementioned two- 
dimensional histogram. This method avoids possible biases 
related to the presence of bright stars near the centre. We de¬ 
rived the cluster ellipticity e running the task ellipse within 
IRAF/STSDA^ on the surface number density images. Fi¬ 
nally, we derived radial surface number densities by dividing 
the area sampled by the images in a series of elliptical an¬ 
nuli, centered on the cluster, and accounting for the spatial 
and photometric completeness in each annulus. 

We only considered stars brighter than the 30% com¬ 
pleteness limit in the core region, corresponding to F438W 
23.5 mag. The outermost data point is derived from the 
AGS parallel observations in a field located ~ 5f5 from the 
cluster centre. The radial surface number density profile was 

® STSDAS is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 
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fitted using a lKin 3 lll962l) model combined with a constant 
background level, described by the following equation: 

n(r) = no ( , ^ - , ^ ^ | + bkg (4) 

\V^ + 0/rc)^ VTT^J 

where no is the central surface number density, is the core 
radius, c = rt/vc is the King concentration index (rt being 
the tidal radius), and r is the geometric mean radius of the 
ellipse (r = a \/l — e, where a is the semi-major axis of the 
ellipse). In Figure [T^ we show the best-fit King model, ob¬ 
tained using a minimization routine. We reported also 
the derived number density values along with other relevant 
parameters. Our derived core radius of = 3.18 (± 0.12) pc 
is significantly larger than the literature value for N GC 1856 
(tc = 1.14 pc. iMcLaughlin fc van der Marell 120051 1. In rec¬ 
onciling this difference, we note that our King-model fit was 
done usi ng completeness-corrected surf a ce num ber densities, 
whereas IMcLaughlin fc van der Marell (l2005l i used surface 
brightness data to derive the structural parameters of the 
cluster. The latter method is sensitive to the presence of 
mass segregation in the sense that mass-segregated clusters 
will appear to have smaller radii when using surface bright¬ 
ness data than when using plain surface number densities. 
To check whether NGC 1856 is actually mass segregated, we 
derived the cumulative completeness-corrected radial distri¬ 
bution of bright and faint stars in the WFC3 image. “Bright” 
stars are selected by the magnitude cut F438W < 20.5 mag, 
whereas “faint” stars are selected in the magnitude range 
21.0 < F438W < 23.5 mag. The obtained cumulative radial 
distributions are shown in Figure 1131 Note that the bright 
stars are clearly more centrally concentrated than the faint 
ones, conhrming that a signihcant degree of mass segrega¬ 
tion is present in the cluster. 

7.2 Dynamical evolution and cluster escape 
velocity 

We estimated the cluster mass and escape velocity as a func¬ 
tion of time, going back to an age of 10 Myr, after the cluster 
has survived the era of violent relaxation and when the most 
massive stars of the hrst generation, proposed to be candi¬ 
date polluters in literature (i.e., FRMS and massive binary 
stars), are expected to start losing significant amounts of 
mass through slow winds. The current mass of NGC 1856 
was determined from its integrated-light F-band magnitude 
listed in Table (2] We determined the aperture correction 
for this magnitude from the best-fit King model derived in 
Sect. Oby calculating the fraction of total cluster light 
encompassed by the measurement aperture. After correct¬ 
ing the integrated-light V magnitude for the missing cluster 
light beyond the measurement aperture, we calculated the 
total cluster mass adopting the values of Av, [Z/H], and age 
listed in Table[2l This was done by interp olation between the 
MJLv values in the SSP mod e l tabl es of iBruzual &: Charlotl 
ll2003lf ■ assuming a ISalpeteJ (119551 1 initial mass function. 
The latter models were recently found to provide the best 
fit (among popular SSP models) to observed integrated- 
light photometry of LMC clusters with ages and metallic- 
ities measured from CMDs and s pectroscopy of ind ividual 
stars in the 0.2-1 Gyr age range (IPessev et al.ll200al . 

We calculated the dynamical evolution of the star clus¬ 


ter following the prescriptions of G14. Briefly, the evolution 
of cluster mass and radius was evaluated with and without 
initial mass segregation, given the fundamental role that the 
latter property plays in terms of the early e volution of the 
cluster’s expansion and mas s loss rate (e.g., iMackev et al.l 
l2008bl : IVesperini et al.ll2009l i. For the case of a model clus¬ 
ter with initial mass segregati on, we adopted the re sults of 
the simulation called SG-Rl in iD Ercole et all (l2008l L which 
involves a tidally limited model cluster that features a level 
of initial mass segregation of re/re,>i = 1.5, where re,>i 
is the effective radius of the cluster for stars with A4 > 1 
MQ(see G14 for a detailed description of the reasons for this 
choice). 

Table [3] lists the derived masses and escape velocities at 
an age of 10 Myr, obtained with and without initial mass seg¬ 
regation. Escape velocities are calculated from the reduced 
gravitational potential 14sc(r, t) = (24>tid(t) — , 

at the core radius. Here 4>tid is the potential at the tidal 
(truncation) radius of the cluster. The choice to calculate 
the escape velocity at the cluster core radius is related to 
the prediction of the “in situ” scenario, where the second- 
generati on stars are formed i n the innermost region of the 
cluster dP’Ercole et al.l lioosll . For convenience, we define 
Ksc, 7 (r) = Ksc {r, t — lO^yr), and refer to it as “early escape 
velocity”. To estimate “plausible” values for Fesc, 7 we used 
a procedure that involves various results from the compila¬ 
tion of Magella nic Cloud star cluster properties and N-body 
simulations by iMackev et al.l ll2008bfl . Briefly, they showed 
that the maximum core radius seen among a large sample of 
Magellanic Cloud star clusters increases approximately lin¬ 
early with log(age) up to an age of ~ 1.5 Gyr, namely from ~ 
2.0 pc at ~ 10 Myr to ~ 5.5 pc at ~ 1.5 Gyr. Conversely, the 
minimumcoie radius is ~ 1.5 pc thro ughout the age r ange 10 
Myr-2 Gyr. Using N-body modeling IMackev et al.l (l2()08bfl 
showed that this behavior is consistent with adiabatic ex¬ 
pansion of the cluster core in clusters with different levels 
of initial mass segregation, in that clusters with the highest 
level of mass segregation experience the strongest core ex¬ 
pansion (see G14 for a more detailed discussion). Under this 
assumption, we derive that the “plausible” early escape ve¬ 
locity for NGC 1856 is 14sc = 17.1±0.8 hms”*^, high enough 
to retain material shed by the slow winds of the polluters 
stars. 

Figure [TT] shows the escape velocity of NGC 1856 as 
a function of age, derived based on the assumed level of 
mass segregation presented above. The critical escape ve¬ 
locity range of 12-15 kms“^ derived by C14 and G14 is 
depicted as the light grey region in Figure 1141 The region 
below 12 kms”*^, representing, as stated above, the velocity 
range in which eMSTOs are not observed in intermediate- 
age star clusters (C14), is shown in dark grey. We note that 
Uesc for NGC 1856 is > 14—15 kms”*^ out to an age of « 100 
Myr. We recall that this is equivalent to the lifetime of stars 
of ~ 5 Mq (i.e., IM-AGB stars, see IVentura fc D’Antonal 
I 2 OO 9 II so that the slow winds of massive binary stars and 
IM-AGB stars of the first generation should have produced 
significant amounts of “polluted” material within that time. 
As mentioned above, this situation is similar to that shared 
among eMSTO clusters of intermediate age in the Magel¬ 
lanic Clouds (G14). It thus seems plausible that a signihcant 
fraction of that material may have been retained within the 
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V 

Aper. 

Aper. corr. 

[Z/H] 


Age 

Tc 

Te 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

10.06 ± 0.15 

31 

0.38 

-0.30 

0.47 

300 

3.18 ± 0.12 

8.00 ± 0.90 


Table 2. Physical properties of the star cluster. Columns (1): Integrated-light V magnitude from iBica et al.l lll996ll . (2): Aperture radius 
in arcsec used for the integrated-light measurement. (3): Aperture correction in mag. (4): Metallicity in dex. (5): Visual extinction in 
magnitude. (6): Mean age in Myr. (7): Core radius in pc. (8): Effective radius in pc. 


Current 

(1) 

log (Mci /Mq) 
10^^ yr w/o m.s. 

(2) 

10^ yr with m.s. 

(3) 

Current 

(4) 

Vesc 

10^^ yr w/o m.s. 

(5) 

(kms 

10^ yr with m.s. 
(6) 

10^ yr “plausible” 

(7) 

5.07 ± 0.07 

5.15 ± 0.10 

5.37 ± 0.10 

11.4 ± 0.7 

13.7 ± 0.7 

20.0 ± 0.7 

17.1 ± 0.8 


Table 3. Dynamical properties of the star cluster. Columns (1): Logarithm of the adopted current cluster mass. (2-3): Logarithm of the 
adopted cluster mass at an age of 10^ yr without(with) the inclusion of initial mass segregation. (4): Current cluster escape velocity at 
the core radius. (5-6): Cluster escape velocity at the core radius at an age of 10^^ yr without(with) the inclusion of initial mass segregation. 
(8) “Plausible” cluster escape velocity at an age of 10^^ yr. 



Age (Myr) 

Figure 14. Escape velocity 14sc as a function of time. The light 
grey region represents the critical range of Fesc mentioned in Sect. 
[3 i.e., 12 — 15 kms”The region below 12 kms”in which I4sc 
is thought to be too low to permit retention of material shed by 
the first stellar generation, is shown in dark grey. ±1 (t errors of 
Fesc are shown by dashed lines. 


potential well of NGC 1856, allowing second-generation stars 
to be formed. 

Finally, we note that the mass and escape velocities of 
NGC 1856 were likely high enough for the cluster to be able 
to accrete a significant amount of “pristine” gas from its 
surroundings in the first several tens of Myr after its birth. 
This would have constituted an additional source of g as to 
form second -generation stars ( [Conroy fc Sperg^l201ll : G14; 
but see also iBastian fc StraclCT 20141 ') . 


8 SUMMARY AND CONCLUSION 

We presented results obtained from a study of new, deep 
HST/WFC3 images of the young (age ~ 300 Myr) massive 


star cluster NGC 1856 in the LMC. After correction for dif¬ 
ferential reddening, we compared the CMD of the cluster 
with Monte Carlo simulations, both of one “best-fit” SSP 
and of two SSPs of different ages, in order to investigate the 
MSTO morphology and to quantify the intrinsic width of the 
MSTO region of the cluster. We studied the physical and dy¬ 
namical properties of the cluster, deriving its radial surface 
number density distribution and determining the evolution 
of cluster mass and escape velocity from an age of 10 Myr 
to 1 Gyr, considering the possible effects of initial mass seg¬ 
regation. 

The main results of the paper can be summarized as 
follows: 

• NGC 1856 shows a broad MSTO region whose width 
can not be explained by photometric uncertainties, LMC 
held star contamination, or differential reddening effects. 
Comparison with Monte Carlo simulations of a SSP shows 
that the observed pseudo-age distribution is signihcantly 
wider than that derived from a single-age simulation. Con¬ 
versely, combining two SSPs with different ages, we obtain a 
set of pseudo-age distributions that reproduce the observed 
one quite well. By considering the luminosity function of 
the RC feature, we conclude that a best ht is achieved with 
a combination of SSPs with ages of 300 Myr and 380 Myr 
and a mass fraction for the younger component of «55%. 
The observed pseudo-age distribution shows two distinct 
peaks with a similar maximum value and an uniform de¬ 
cline towards younger and older ages, with respect to the 
peaks. However, the small separation in age between the two 
peaks prevents us to conclude whether the morphology of 
the MSTO can be better explained with a smooth spread in 
age or by two discrete bursts of star fo rmation. These results 
do no t agree with the conclusions of IBastian fc Silva-Villal 
ll2013t) who conclude that the CMD of NGC 1856 is consis¬ 
tent with a single age to within 35 Myr. We expect that this 
difference is due to the fact that our data are signihcantly 
deeper, name ly ^ 6 mag. Consequently, w e suggest that the 
arguments of IBastian fc Silva-Villal (l2013f) against the “age 
spread” scenario should be considered with caution. 
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• We use V, I, and Ha images to select and study candi¬ 
date (“putative”) pre-MS stars in NGC 1856. The numbers 
of “putative” pre-MS stars in the cluster field and the con¬ 
trol (background) field are found to be similar, suggesting 
that the detections can be considered spurious and/or asso¬ 
ciated to residual noise. This indicates that star formation 
is not currently ongoing in the cluster. 

• The “stellar rotation” scenario for the nature of the 
eMSTO phenomenon has been tested by evaluating rotat¬ 
ing and non-rotating isochrones of the same age with the 
DRC CMD. Overall, a reasonable range of rotation velocities 
seems to be able to reproduce the MSTO properties quite 
well, albeit only if the rotational mixing efficiency is signifi¬ 
cantly higher than typically assumed values. However, several 
properties of MSTOs and RCs among eMSTO star clusters 
in the age range of 1 - 2 Gyr are inconsistent with such high 
rotational mixing efficiencies. This seems to indicate that 
a range of stellar ages provides a more likely explanation 
of the eMSTO phenomenon (including the wide MSTO of 
NGC 1856) than does a range of stellar rotation velocities, 
altough, with the current data alone, the latter hypothesis 
cannot formally be discarded. In this context, new stellar 
rotation measurements in young and intermediate-age star 
clusters, combined with new models of theoretical tracks and 
isochrones for rotating stars will be of fundamental impor¬ 
tance to address the role of stellar rotation in the eMSTO 
phenomenon. 

• The dynamical properties of NGC 1856 derived from 
our data suggest that the cluster has an early escape veloc¬ 
ity of ~ 17 kms“^, high enough to permit the retention of 
material shed by the slow winds of polluter stars (IM-AGB 
stars and massive binary stars). The cluster escape velocity 
remains above the threshold value of 14-15 kms“^ for « 
80 -100 Myr, long enough for the slow winds of IM-AGB 
stars to have ejected their envelopes. This material would 
likely have been available for second-generation star forma¬ 
tion which could have caused the wide MSTO. Furthermore, 
these early escape velocities are consistent with observed 
wind velocities of ejecta of IM-AGB stars and with those 
seen in massive binary star systems. 
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APPENDIX A: PHOTOMETRIC CATALOG 

In Table [Q we report positions and magnitudes, with the 
associated errors, for the first ten objects in our final pho¬ 
tometric catalog. The complete table of stellar photometry 
is available upon request. 
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ID 

(1) 

X 

(2) 

Y 

(3) 

F438W 

(4) 

err 

(5) 

F555W 

(6) 

err 

(7) 

F656N 

(8) 

err 

(9) 

F814W 

(10) 

err 

(11) 

1 

1277.992 

16.072 

27.066 

0.280 

99 

99 

99 

99 

23.549 

0.195 

2 

1279.163 

19.076 

99 

99 

27.625 

0.262 

22.958 

0.317 

99 

99 

3 

1274.036 

20.648 

26.536 

0.156 

99 

99 

22.620 

0.292 

99 

99 

4 

1300.449 

22.857 

27.221 

0.289 

99 

99 

22.664 

0.239 

99 

99 

5 

1319.126 

26.275 

21.564 

0.019 

99 

99 

20.680 

0.072 

99 

99 

6 

1296.704 

26.879 

27.217 

0.377 

99 

99 

22.607 

0.315 

99 

99 

7 

1306.224 

27.537 

27.096 

0.208 

99 

99 

22.816 

0.218 

99 

99 

8 

1269.086 

28.177 

27.537 

0.248 

99 

99 

23.518 

0.462 

99 

99 

9 

1290.509 

29.335 

25.467 

0.076 

24.877 

0.058 

99 

99 

23.678 

0.041 

10 

1278.772 

31.292 

24.668 

0.059 

24.051 

0.045 

99 

99 

23.047 

0.034 


Table Al. Stellar photometry. (1): Identification number. (2-3): X and Y coordinates in pixels. (4-11): magnitudes and associated errors. 







